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Abstract 

We present numerical simulations of colliding wave packets in spontaneously 
broken SU(2) Yang-Mills-Higgs theory. Compared with pure Yang-Mills the- 
ory, introducing the Higgs field leads to new aspects in the dynamics of the 
system. The evolution of the gauge field and the Higgs field is investigated 
as a function of the amplitude of the wave packets and of the mass ratio of 
the Higgs and the gauge boson. We find regions in our parameter space in 
which initial wave packets scatter into final configurations with dramatically 

different momentum distributions. 
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I. INTRODUCTION 



Collisions between classical wave packets have recently been studied numerically for sev- 
eral interacting relativistic field theories |l]-f|] . Interest in this topic arose in connection with 
expectations that the rate of multiparticle production processes in electroweak interactions, 
which can manifest themselves, e.g., in baryon number violation, might be unsuppressed at 
high energies |J. 

The non-perturbative nature of the baryon number violating amplitude |3| demands a 
corresponding non-perturbative approach as provided by semiclassical techniques. The main 
difficulty in semi-classical approaches is the treatment of the 2 — > many particle transition 
amplitude, since the initial state of high energy particles is not semiclassical at all, and loop 
contributions are essential in general (see, e.g. 0). Possible techniques for circumventing 
this difficulty of the semiclassical approach have been proposed and studied in the literature 




In the classical approach to scattering, the question is: does there exist a mechanism for 
energy transfer from high frequency modes, which corresponds to two high energy particles, 
to low frequency modes representing a multiparticle final state? At first glance, the answer 
to the question, formulated in terms of nonlinear dynamics, seems to be affirmative since the 
gauge field theories are nonlinear. However, the studies of (1 + l)-dimensional abelian Higgs 
model and A<^ 4 -theory j2[ have shown no indication for a non-perturbative mechanism 
providing the coupling between the initial high and the final low frequency modes. For ex- 
ample, in [p]] the wave packets always passed through each other without being destroyed. It 
is important to note that in || the initial states were always chosen to have small amplitude, 
which made the nonlinear terms less important. 

The important issue here is that the results are strongly influenced by the nonlinearity 
due to the non-abelian spin-field coupling, which is absent in abelian models. It is this 
coupling that is responsible for the infrared instabilities of the pure non-abelian gauge theory. 
This can be seen from the linearized equation describing small perturbations a c around an 
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SU(2) background field (in background gauge): 

{EfaY - 2g e abc F b ^ = , (1) 

where = — ig[A^, ] is the gauge covariant derivative. The second term in (|IJ) may 
have any sign. In particular, this essentially non-abelian coupling drives an instability for 
perturbations with isospin polarization orthogonal to the isospin of a standing wave, which 
leads to a growth of low frequency modes from initial high frequency modes |J. This may 
imply the existence of classical trajectories of the type required for multiparticle production 
0, if this instability persists in more realistic case, e.g., collision of localized gauge field 
wave packets. 

From a more general point of view, the observed inability of the nonlinearity to furnish 
a mechanism for the formation of strongly inelastic final states is, in our opinion, intimately 
connected with the integrable nature of the classical systems considered in |I[ and [Q. It 
is well-known that non-abelian gauge theories are non-integrable in the classical limit , and 
exhibit dynamical chaos |TT| , P^| Q This dynamical stochasticity of the non-abelian gauge 



fields together with their mentioned dynamical instability are possible sources of the non- 
perturbative mechanism for the coupling between high frequency modes and low frequency 
modes. At the same time, it is important to recall the special role of the Higgs field as a 
mechanism for the suppression of the dynamical chaos of the non-abelian gauge fields [0 . 

With this in mind, we studied f| the collision of two SU{2) gauge field wave packets, 
homogeneous in the transverse plane. As we expected, based on previous results ||, the 
collision of essentially non-abelian initial configurations trigger the decay of initial high 
frequency modes into many low frequency modes with dramatically different momentum 



1 Strictly speaking, chaos only sustains for solutions of finite energy density [12|. Finite energy 
solutions in (3 + l)-dimensions will spread out in space at late times and will therefore linearize. 
However, nummerical results indicate that at intermediate times these fields generally exhibit 
exponentially growing perturbations. 



distributions, whereas for abelian configurations (parallel relative isospin polarizations) wave 
packets pass through each other without interaction. 

The present paper, which studies collisions of wave packets in the 577(2) Higgs model, 
is a generalization of the earlier work in two directions. Firstly, it is an extension of to 
the non-abelian Higgs model. Secondly, it is a generalization of the previous work [|J to 
the case where the SU(2) gauge symmetry is spontaneously broken by a Higgs field and the 
fundamental excitations of the gauge field are massive. 

One expects that the explicit mass scale introduced by the Higgs field will act as a cut- 
off on the low frequency excitations, potentially leading to drastic changes in the coupling 
between high and low frequency modes. We will see that the real situation is more compli- 
cated, and the ratio \/g 2 of the Higgs self-coupling A and gauge coupling g and the vacuum 
expectation value v of the Higgs field are essential parameters. 

This paper is organized as follows. In Section II, we formulate the problem. In Section 
III, we present results from our numerical simulations and discuss their implications. Section 
IV is devoted to an extended discussion of our results. In Section V, we conclude and indicate 
possible directions for future research. 



II. FORMULATION OF THE PROBLEM 

In this section, we describe scattering of classical wave packets in the non-abelian Higgs 
model and its numerical formulation on the lattice. In particular, we give a brief discussion 
of the scaling properties of the classical dynamics. This section is based on previous work 



A. The non-abelian Higgs model 

Here we give a brief discussion on the spontaneously broken SU(2) Yang-Mills theory, in 
which a charged scalar isodoublet field, the Higgs field, is coupled to the gauge field. This 
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model retains the most relevant ingredients of the electroweak theory. The action describing 
this model in (3 + l)-dimensions is given by 



S = J d?xdt -\r(F^F» u ) + Ur ((D^D^) - A Qtr($ t $) 



(2) 



with Dp = — igA a r a /2, F, u = F a r a /2 = (i/g) [£>„, D v \ and 



$ = 0° - i T a (j) a 



(3) 



where r a (a = 1,2,3) are Pauli matrices. Following the notation of fl4 l, we represent the 
complex Higgs doublet by a quaternion, which is convenient for numerical manipulation. 
Clearly, this maintains the correct number of degrees of freedom in the Higgs field. 
By a scaling transformation 



^ = gvx„ $' = ®/v, A'^ = A„/v , 



(4) 



we obtain the action in terms of the primed quantities 

S = (1/g 2 ) f d 3 x'dt' --tr^F'^) + ~tr ((D^')^'^') " A ' f-tr($' t $ / ) - 1 
J 2 2 \ 2 



(5) 



with A' = \/g 2 . 

Within classical physics, the prefactor 1/g 2 in (|5|) is irrelevant, leaving A' = X/g 2 as the 
only relevant parameter in the action. Note that A' is proportional to the square of Mff/Myy, 
the mass ratio of the Higgs and H^-boson. 

The elementary excitation modes p and are best described in the unitary gauge 



$ = (v + p/V2)U(6) , 
A, = U{9)W,U-\e) - (1/ig) (d,U(6)) U-\6) 



(6) 
(7) 



where U{6) = exp(ir a 9 a ), p describes oscillation of the Higgs field about its vacuum expec- 
tation value, and is the field of the gauge boson. and p obey the following classical 
equations of motion 
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[£>„, F^\ + M 2 W W V + ±=g 2 vpW v + \g 2 p 2 W» = , (8) 
(«9 M <9 M + M|)p + 3V2Xvp 2 + Ap 3 - ^g 2 vW;W a ^ - ^g 2 pW^W afl = , (9) 

where M H = 2v\f\ and Mvk = W^2. ^ and are denned in terms of W^. After a 
scaling transformation similar to that in (f|), it is easy to see that the above equations of 
motion depend on a single parameter: the mass ratio Mh/Mw- However, in the simulation 
of colliding wave packets, there are other parameters involved in the initial condition. 



B. Scattering of classical wave packets 

Our numerical study is based on the Hamiltonian formulation of lattice 577(2) gauge 
theory [HJ (see [fi|,|T6| for more details), in which the dynamic variables are link variables 
defined as 

U t = exp(-igaA c e r c /2) , (10) 

where I stands for the link index. As in we work on a one dimensional lattice with a 
physical size L = Na, where iV is the number of lattice sizes and a the lattice spacing. We 
arrange initially two Gaussian wave packets with average momenta k = (0, 0, k), and width 
Ak. Our goal is to simulate the collision of two H^-boson wave packets in the background 
of the Higgs condensate. 

Before actually constructing the wave packets, one has to deal with the gauge fixing 
problem. In the Hamiltonian formulation of lattice gauge theory, temporal gauge Aq = 
is most convenient. On the other hand, one must construct the jy-boson wave packets in 
the unitary gauge and then transform back to obtain the initial conditions in the temporal 
gauge. The gauge field for a configuration of two well-separated wave packets in the unitary 
gauge can be written as 

W c '» = + , (11) 
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with c being the isospin index and /i the Lorentz index. is a left-moving wave packet, 

initially centered at Zl, is a right-moving wave packet, initially centered at zr. In our 
simulation, zi and zr are chosen in a way such that the two wave packets are positioned 
symmetrically about the center of the lattice. Specifically, we take transversely polarized 
wave packets 

W^ = (0,0,l,0)n c L ^(z-z L ,-t) , (12) 
W c jt = (0, 0, 1, 0)n c R ^(z - z R , +t) , (13) 

with n c L and n c R being the polarization vectors in isospin space. We choose n c R = (0, 0, 1) 
fixed and leave n c L free to be varied. The above choice satisfies the relation d^W^ = 0. 

Because we have chosen transversely polarized wave packets^, they already satisfy the 
temporal gauge condition Wo = 0. If we had chosen, instead, longitudinally polarized wave 
packets (as in ref. |l|), we would have needed to apply a gauge transformation ?7(x, t) which 
transforms to 

A, = UW.U- 1 - (l/tg)(d ll U)U- 1 , (14) 

such that A = 0, i.e., 

d U = igUW . (15) 

In our case, obviously, U — 1. 

To construct the wave packets, we need to specify the functional form of ^(x, t). A 
right-moving wave packet centered at z = at time t = with mean wave number k, width 



Ak, and mean frequency u = Jk 2 + is described by 



2 Note that under realistic conditions the luminosity for transversely polarized gauge bosons in 
proton-proton system is typically two orders of magnitude higher than for longitudinally polarized 
ones and increases with energy Jl7| . 
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iP(z,t) 



dk e -^-w/mk? 



e i(u(k z )t-k z z) _|_ c c 



(16) 



with uj(k z ) = Jk z + Myy and 



Q = uj 



1 + 
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1 



1 

UJ 



'AkJ + f;Ak^ 



UJ 



(17) 



where the amplitude of the wave packet is fixed by requiring energy equal to hu ("one 
particle") per cross sectional area a. In the following we will set % — 1. 
Performing the /^-integral at t — gives 



2Ak 



t=o 



-(Afc,z) 2 /2 



cos(/c,2) 



Since the differential equations are of second order in time, one also needs to specify ip 
dijj/dt at t = 0, which is found to be 



(18) 



■<P\t=o 



AiryJnVLAka 



dk z uj(k z )e 



-(fc z -fc) 2 /2(Afc) 2 



e i(ui(k z )t~k z z) _ c c 



t=0 



2toAk 



yJlia{Vt/uj) 



1 

+ 2 



-(Afc 2 ) 2 /2 



P 2 ' 



sinf&z) + — ^— A/cz cos(fc;z) 
— J [1 - (Ah) 2 ] sin(A;z) + O (( — 



(19) 



Furthermore, the initial condition for the Higgs field is chosen as the vacuum solution: 



iP = v, & a = 0, dp = <b a = at t = 



(20) 



To determine the number of independent parameters, we make use of the scaling trans- 
formation (H) for ip\t=o- In terms of the primed quantities, equation (fL8|) reads 



V>'|tf=o = (l/f)^k=o 
= (V- 1/4 



Ak/M] 



w 



\ (n/M w )(aM^/g^) 



-(Ak/M w )*z>*/i cos [(fc/M Ty )// V / 2], 



where 



Mi 



1 + (fc/M^) 2 



1 + - 



1 (Afc/Mi 



4 [1 + (fc/M^) 2 ] 2 



+ C?((Afc/fc)' 



(21) 



(22) 
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The above initial condition contains three dimensionless parameters: k/Mw, Ak/Mw, and 
aM w /g 2 . There appears one more parameter in the initial condition, i.e., the relative 
rotation in isospin space between the two wave packets, which we denote by 6 C . 

Combining equations of motion and initial condition, our ansatz has five independent 
parameters: M H /M W , k/Mw, Ak/M w , aM w /g 2 , and 9 C . The parameter k/Mw, which sets 
the energy of collisions in units of the H^-boson mass, is referred to as the energy parameter. 
Ak/Mw, whose inverse specifies the width of each wave packet in position space, can be 
called the width parameter. The amplitude of each wave packet depends on k/Mw as well 
as Ak/Mw, but more crucially, on aM w /g 2 . In our simulation, we always require k S> Ak 
so that the wave packets are well-defined objects. Furthermore, we choose k 3> M w to 
model high energy scattering. 

In our numerical calculations, the SU(2) coupling constant g was fixed to be 0.65. How- 
ever, due to the scaling properties of the equations of motion and the initial conditions, the 
results of our calculation do not depend on the particular choice of g and v. This can be 
verified by the fact that the amplitude of the wave packets only depends on the ratio of a 
and g 2 . Also, since the dynamics does not depend on a particular choice of a or g as long 
as the ratio u / g 2 is fixed (assuming that Mw and other parameters remain fixed), we can 
predict, from the result for one coupling g at a certain value of a, the result for another 
coupling g' at a different a' = (g' /g) 2 a. Hence, a change of the value for the gauge coupling 
g simply corresponds to a rescaling of the parameter a controlling the amplitude of initial 
wave packets. 

III. NUMERICAL RESULTS 

A. Dependence on the mass ratio Mh/Mw 

As established in the previous work Q, the behavior of the wave packet collisions is 
governed by the nonlinearity due to the self-interaction of the gauge field. For two wave 
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packets of parallel isospin polarizations in the pure Yang-Mills theory, where the nonlinear 
self-coupling in the gauge field is absent, we found no indication of any interaction [Q. This 
provided a check on our numerical procedure and showed that the artificial interactions 
introduced by the formulation in terms of compact lattice gauge fields did not affect the 
results. 

Here, for the Yang-Mills-Higgs system, the situation is more involved. Besides the non- 
linearities due to the gauge field self-interaction, there exist other nonlinearities induced by 
the gauge-Higgs coupling and by the Higgs self-interaction. However, in the case where most 
energy remains contained in the gauge field and the Higgs field is only slightly excited, one 
can expect that the gauge field self-interaction will be the major contributor to the non- 
linearities observed in the system. Note that the amplitudes of the gauge and Higgs fields 
shown in the figures below are in the temporal gauge. This means that the longitudinal part 
of the gauge field is not fully represented in the figures. 

The top rows of Figure 1 and Figure 2 show a few "snapshots" of the space-time develop- 
ment of the colliding VT-boson wave packets with Mjj = Myy = 0.126, a = 0.336, k = 7r/5, 
and Ak = 7r/100 for parallel (Figure 1) and orthogonal (Figure 2) isospin polarization, re- 
spectively. The figures show the absolute magnitude of the scaled gauge field amplitude, 
| A' | = \A\/v. For parallel isospin orientations, the result of the "collision" is a slight distor- 
tion of the initial wave packets showing no sign of significant inelasticity. In contrast, the 
top row of Figure 2 illustrates that the collision of two wave packets with orthogonal relative 
polarizations in isospin space is strongly inelastic. 

The difference between the two figures is even more clearly illustrated by looking at the 
evolution of the absolute value of the Fourier transform of the gauge invariant energy density 
(scaled by v 2 ) 

£{k,t) 1 



£(M) 



Jd 3 x e ikx Tr [E 2 (x,t) + B 2 (x,t)] , (23) 



v 2 Av 2 

where E is the gauge electric field and B the gauge magnetic field. It is seen that the 
spectrum for the parallel isospins (median row in Figure 1) does not change its shape dra- 
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matically, while for the case of orthogonal isospins (median row in Figure 2), the spectrum 
spreads out widely. The spike at k = in these spectra corresponds to the total energy 
contained in the transverse gauge field. From its slight decrease in Figure 2, we see that the 
energy transferred to the Higgs and the longitudinal gauge fields during the collision is small 
(~ 10%). This reflects the fact that the nonlinearity due to the gauge field self-coupling 
dominates. Furthermore, the bottom rows of Figures 1 and 2 illustrate the time evolution 
of the Higgs field excitations^ around its condensate value v (which is scaled to unity) ac- 
companying the collision process shown in the top rows of Figures 1 and 2. Here we have 
plotted the square |$'| 2 = |$| 2 /t> 2 of the Higgs field as a function of space coordinate at 
three different times. Throughout our simulations, we have kept k min <C Ak <c; A; <C k max , 
where k m i n and k max are the minimum and maximum momentum on the lattice, respectively. 
This ensures that the wave packets are smooth on the lattice. But during collisions, unlike 
quantum mechanics, classical dynamics does not provide a mechanism for stopping power 
from flowing to very low frequency modes (close to k min ) or to very high frequency modes 
(close to k max ). In our calculations, the power flowing to high frequency modes does not 
cause a deterioration in the local smoothness of the gauge field at the end of the simulations. 

To display dependence on the mass ratio r = Mn/M\y, the top row of Figure 3 shows the 
collision of two orthogonally polarized W wave packets at the final time (t = 580) for three 
different values of r = MujM^/. It is seen that the "inelasticity" is more pronounced for 
small r. On the other hand, the distortions in the wave packets still survive at large r (even 
at r = 100, not shown here). This can be understood as follows. Remember that there are 
three sources of nonlinearity, namely, gauge field self-coupling, gauge-Higgs coupling, and 
Higgs self-coupling. As the Higgs mass increases, the Higgs modes begin to decouple. As 
a result, the interaction between gauge and Higgs field diminishes and hence contributes 



3 As seen from (||) oscillations of the gauge boson field act as a source for Higgs field excitations. 
The equation (|8|) for the gauge field does not possess a source term. 
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less to the nonlinear effects. The gauge field self-interaction is not affected by the change 
in the Higgs mass and acts as the main contributor of nonlinear effects observed during 
the collision. The median row of Figure 3 is the analogue of the top row of Figure 3 in 
momentum space, as defined in (pci|). Again, it gives a clearer picture of the inelasticity of 
the collision process. The bottom row of Figure 3 demonstrates that the amplitude of the 
Higgs field excitations becomes smaller as the Higgs mass is increased, while their frequency 
increases with Mh- 

It is remarkable that for small Higgs mass, as seen for r = 0.1 in Figure 3, the Higgs field 
oscillates not about the vacuum expectation value v but rather about zero. The observed 
behavior holds even at larger values of r, up to r « 0.5 (not shown here). This suggests 
that for not too large r, the collision of gauge field wave packets, accompanied by energy 
transfer from gauge field to Higgs field, leads to restoration of the broken symmetry. This 
phenomenon occurs for gauge field configurations with large amplitude (]. Indeed, it is easy 
to see that p = —\[2v (i.e. |$| = 0) is an exact solution to equation (|9|). Inserting p = —\/2v 
into equation (|§p leads to the pure Yang-Mills equation for massless VT-bosons. In terms of 
excitations around this state, x — P + V%v, we rewrite equations (§) and (|]) as 

[D li ,F^ + ±g 2 X 2 W v = 0, (24) 



d^-M 2 H {\ 



g 2 W 2 



X + Ax 3 = , (25) 



l\v 2 

where W 2 = — (W"") 2 < for transverse polarized wave packets (the sum over spatial index 
i and isospin index a is assumed here and in below). After dropping a constant term At> 4 , 
the corresponding effective potential describing the excitations x nas the following form: 

V( X , Wg = -\v 2 (l - V ) X 2 + ^X A , (26) 



^The idea that the restoration of vacuum symmetry is possible in the background of intense gauge 



fields was first noted in [18]. Also see [|ij|]), where the role of external gauge fields in the restoration 



of broken symmetries was considered. 
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where we denote 

_ g 2 (W?) 2 1 (W? 



]\v 2 r 2 

ra\2 



(27) 



In the following we use 77 as a parameter in which the true intensity (Wf) of the high 
frequency gauge field pulses is replaced by its space-time average (W 2 ). Depending on 
whether r\ < 1 or 77 > 1, the potential (GBI) has two different stable minima: 



for 77 < 1, xmin = ±V2v{l - V ) 1/2 , i.e. |$| = v{l - r/) 1/2 , (28) 
for 77 > 1, Xmm = , i.e. |$| = . (29) 

Stable excitations about these "vacua" have the following squared masses: 

M% = M^(l - 77)0(1 - r,) (30) 

~ M 2 

M 2 H = ^(1- V )[1 + 6(1- V )] (31) 

Thus for 7] > 1 , the broken symmetry is restored and oscillations of the scalar field occur 
about the symmetrical state |$| = 0, not about |$| = v. The effective mass of the gauge 
bosons in the region where the symmetry is restored vanishes. For 77 < 1, the ratio between 
the effective masses Mh and My/ has no dependence on 77 and remains r = Mh/Mw. 
Relations ( |3"0"| ) and ([IT]) are characteristics of a second order phase transition. The expression 



(p7|) for 77 shows that in the regime of large (W? - ) 2 > v 2 , this phase transition can occur for 
experimentally favorable mass ratio r > 1. 

Oscillations of the scalar field around the new symmetrical minimum |$| = are clearly 
seen for r = 0.1 where 77 m 22.3 (see Figure 3). These numerical results provide indications 
for transition from the phase with spontaneously broken SU(2) symmetry and asymmetric 
vacuum to the phase with restored SU(2) symmetry and symmetric vacuum as a result of 
the collisions. Since rj is the only relevant parameter in question here, this transition can 
occur either for small A (light Higgs) or for large A (heavy Higgs) if amplitude of the gauge 
field is large enough. In Figure 3, it is also interesting to notice that the spatial region 
showing symmetry restoration seems to be wider than the region where the colliding wave 
packets stay visibly large. 
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It is clear that the phenomenon discussed above does not depend on the one- 
dimensionality of space in our calculations. However, it is to be expected that the symmetry 
restoration would not persist as long in three dimensions as the wave packets disperse more 
rapidly after the collision causing the squared amplitude (W 2 ) to decay more rapidly. 



B. Yang-Mills and BPS limit 

In the light of the previous work ||, it is instructive to study the limiting case of the 
present system as Mh, My/ — > while fixing r = Mh/Mw (which is chosen to be 1 here). 
Of course, this corresponds to the limit v — > 0, where one expects that the gauge field in 
the Yang-Mills-Higgs system behaves most closely to that in a pure Yang-Mills system. In 
Figure 4, the top row shows snapshots of the collision of two orthogonally polarized wave 
packets with Mh = My/ = 0.001 and the bottom row exhibits the corresponding spectra. 
Qualitatively, these figures show that the time evolution is very similar to that seen in the 
pure Yang-Mills system (see Figures 2 and 4 in [|]]). 

The second interesting limit is the Bogomolny-Prasad-Sommerfield (BPS) limit PH 



where Mh — > but My/ is finite (A — > 0, v fixed). Figure 5 shows the snapshots and 
spectra for this case with My/ = 0.126 and r = 0.01. Again, we display only the orthogonal 
case, which reveals complete destruction of the wave packets as in the pure Yang-Mills limit. 
It is interesting to note that in this limit the static force between equally charged iy-bosons 
vanishes due to the precise cancellation between the photon and (massless) Higgs exchange 
diagrams, which is a result of the electromagnetic duality ||21|| . For a pair of oppositely 



charged W^-bosons, the contributions from these diagrams add to each other, doubling the 
attraction. 



C. Dependence on the initial amplitude and energy 

In our simulation, the most crucial role is played by the initial amplitude of the wave 
packets. As in the pure Yang-Mills case [f§] , we find that the amount of "inelasticity" observed 
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in the present system is closely linked to the magnitude of the dimensionless amplitude 
(pl|). This amplitude depends on several independent parameters: aM^/g 2 , Ak/M w , and 
fl/Mw ~ k/Mw (for k ^> Mw), each of which has a different physical meaning. Noting 
that the gauge coupling constant is fixed to be g = 0.65 throughout our simulation and 
M w is a fixed quantity in reality, the best way to study the amplitude dependence is to 
vary a without changing anything else. By varying a, we find that the nonlinear effects 
increase with amplitude. For a very large a (very small amplitude), we find no indication 
of "inelasticity" in the colliding wave packets at the end of the simulation. 

To search for the energy dependence, we have to study the dependence on k, which 
determines the energy of the initial wave packet. In the meantime, we fix a, My/, and 
Mn. My/ and Mu are chosen to be much smaller than k to model high energy scattering. 
Furthermore, as we change k, Ak is either fixed or changed proportionally to fix the ratio 
Ak/k. In Figure 6, we display snapshots of collisions at the final time for different sets of k 
and Ak. Figure 7 shows the corresponding spectra. The orthogonal isospin cases in Figures 
6 and 7 are shown in the left column with their parallel isospin counterparts in the right 
column. 

Clearly, the observed nonlinear effects in the orthogonal cases are qualitatively similar 
for different k or Ak; while in the parallel cases, regardless of Ak, the nonlinear effects 
disappear as k is increased from 7r/25 to 7r/5. 



IV. DISCUSSIONS 



A. Amplitude Dependence 



Our numerical calculations show for a wide range of parameters that the wave packet 
collisions with orthogonal isospin orientation are strongly inelastic if the scaled initial am- 



plitude (see equation (0), k ^> M w ) 

( 2 Ah y /2 

\ y/nkav 2 J 
15 



(32) 



is of the order of unity. We recall that the expression (|2T|) for the scaled amplitude was 
determined by the condition that the wave packet contains one particle per transverse area 
a. 

Although, in the strict sense, our configurations describe wave packets which are infinitely 
extended in the transverse direction and hence contain infinitely many particles, only a 
finite transverse area influences the dynamics over a finite period of time. As argued in [|J , 
causality restricts that area to cr(T s ) = 7rT s 2 where T s is the elapsed time after the impact of 
the two wave packets. The relevant number of particles in the initial state is therefore given 
by 

= ajn = .Jf 
a a 

What is the lower bound on Nf under realistic conditions, for which strongly inelastic 
events occur? Let us first estimate the constraints on our parameters from a realistic point 
of view. Clearly, we must have k ^> v. The natural spread of any P^-boson wave packet 
produced in high energy interactions is of order Ak ~ v in the comoving reference frame. 
Therefore the typical transverse area of the IV-boson wave packet is of order a fa 1/v 2 . Due 
to Lorentz contraction, its longitudinal momentum spread will generally be much larger than 
v, or of order Ak\\ ~ jv , where 7 fa k/My/ is the Lorentz factor. As a result, Ak\\/k will be 
approximately independent of the collision energy, with a value not much smaller than one. 



Assuming, e.g., Ak/k fa 0.5 in (32), we obtain an amplitude of order unity implying that a 
few VT-bosons per area a in the initial state could produce strong inelasticity. Of course, the 
precise lower bound on the particle number will depend on the detailed shape of the wave 
packets and requires a full three-dimensional analysis. But our estimate shows that strongly 
inelastic events are not excluded for collisions of wave packets containing few particles. In 



this respect, the results of our analysis correspond to those of Singleton and Rebbi JTU| who 
found that few-particle initial states may not be excluded for baryon number nonconserving 
processes resulting in multi-particle final states. 

As mentioned above, the finite transverse size of order v _1 limits the applicability of our 
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calculation for the real three-dimensional case to times T s < v^ 1 . Since the inelasticity is 
clearly revealed for times T s ~ 100 (see e.g. Figure 5), our numerical results apply most 
confidently to very high energy where k/v > 10 2 . 

B. Energy Dependence 

We now turn to the question of the energy dependence of the nonlinear effects seen in the 
wave packet collisions. In the (1 + l)-dimensional abelian Higgs model, the nonlinearities 
were clearly found to decrease with energy [l|]. For the non-abelian Higgs model discussed 
here, the answer is given in Figures 6 and 7. The inelasticity seen in the orthogonal isospin 
cases does not change significantly with energy, while it dies out in the parallel isospin cases 
as k increases. This shows the fundamental role of the non-abelian nature of the W-W 
interaction in the formation of strongly inelastic final states. 

From Figures 6 and 7 one can also see that the inclusion of the Higgs field produces 
new phenomena which are not seen in the pure Yang- Mills system: For initial configurations 
with parallel isospin, in which case nonlinear interactions of the gauge bosons are absent, 
lowering of the parameter k leads to inelastic final states. This is exclusively due to the 
Higgs field since the collision of wave packets with parallel isospin configurations in the pure 
Yang-Mills theory always leads to elastic final states independently of k Q . Of course, this 
pattern indicates that the influence of the Higgs coupling to the gauge field increases for 
smaller k. How does one understand this behavior? For this purpose, we recall that all 
our calculations are in the regime of high energy (k 3> v, M\y)- For the highest energy of 
the collisions (k = 7r/5) where parallel polarized wave packets scatter elastically, one may 
think that the transversely polarized W-W scattering (elastic in these collisions) proceeds 
via exchanges of the gauge and of the Higgs bosons in the tree approximation. The first 
contribution prevails at high energy, but it does not contribute to the inelastic final states for 
the parallel isospin orientation. Inelasticity may arise here only from the nonlinear coupling 
of the gauge and Higgs fields or, in diagrammatic language, due to the Higgs exchange whose 
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contribution increases with the lowering of k. 



V. CONCLUSION 

We have numerically studied collisions between classical wave packets of transversely 
polarized gauge bosons in the spontaneously broken Yang-Mills-Higgs theory. Our main 
results are the following: 

1. We have found evidence for the creation of final states with dramatically different 
momentum distributions (strongly "inelastic" events) for a wide range of the essential 
parameters. 

2. These inelastic events persist at the highest investigated energies (k/M w ~ 10 2 ) for 
collisions with orthogonal isospin polarization, reflecting the essentially non-abelian 
character of the interaction. For parallel isospin configurations, in contrast, the inelas- 
tic events, which are solely due to the Higgs field, occur only for lower energies. 

3. Under more realistic conditions as discussed in Section IV.A, the inelastic events are 
not excluded for initial configurations with few particles. 

4. We have observed, at least for r < 0.5 (with fixed amplitude for the gauge field wave 
packets), the phenomenon of symmetry restoration as a result of the wave packet 
collisions. This transition from the asymmetrical state to the symmetrical one is 
governed by a single parameter rj, which depends on both the mass ratio r and the 
amplitude of the wave packets. 

Summarizing, we conclude that the introduction of the Higgs field (in the broken sym- 
metry phase) does not in general spoil the inelasticity of the final state in collisions with 
orthogonal isospin orientation. Our results provide a strong motivation for exploring related 
phenomena in (3 + 1) dimensions. This would allow one to study collision between real- 
istically shaped wave packets and investigate the particle number content of inelastic final 
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states. Last, but not least, it would be interesting to study the winding number change 
associated with these collisions. 
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FIGURES 

FIG. 1. Collision of two W wave packets with parallel isospin polarizations. We choose 
M H = M w = 0.126, k = vr/5, Ak = vr/100, g = 0.65, and a = 0.336. This simulation, as 
well as all others below, was performed on a lattice of length L = 2048 and lattice spacing a = 1. 
The top row shows the space-time evolution of the scaled gauge field amplitude |A'|, the median 
row exhibits the corresponding Fourier spectra of the gauge field energy density, eq. (|23|), and the 
bottom row shows the space-time evolution of the scaled Higgs field |<F| 2 . The abscissae of top 
and bottom rows are labelled in units of the lattice spacing, and the abscissa of the median row is 
in units of 7r/1024. 

FIG. 2. Same as Figure 1, but for orthogonal isospin polarizations. 

FIG. 3. Mass ratio r = Mh/Mw dependence of the collisions shown for three different values 
of r at the end of our calculation (t = 580). Here we have used orthogonally polarized W wave 
packets. Except the mass of the Higgs, all the other parameters are fixed to be the same as in 
Figure 1. 

FIG. 4. Collision of two orthogonally polarized wave packets in the Yang-Mills limit. Except 
for Mh = Mw = 0.001, all the other parameters are the same as in Figure 1. The top row shows 
the space-time evolution of the scaled gauge field amplitude. The bottom row shows the Fourier 
spectra of gauge field energy density, as defined in (pffi). 

FIG. 5. Collision of two orthogonally polarized wave packets in the BPS limit. Here we fix 
Mw = 0.126 but choose a small mass for the Higgs: Mjj = lO~ 2 M\y- All the other parameters 
are the same as in Figure 1. The top row shows the space-time evolution of the scaled gauge field 
amplitude. The bottom row shows the Fourier spectra of gauge field energy density, as defined in 
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FIG. 6. Final states (t = 1100) of the scaled gauge field for three different sets of k and Ak are 
shown. Here a = 0.504, Mh = Myy = 0.01, and all the other parameters are the same as in Figure 
1. Left column: orthogonal isospin orientations. Right column: parallel isospin orientations. 

FIG. 7. Energy spectra for the cases shown in Figure 6. 
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